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A  recent  RAND  report,  Close-In  Fallout, ^  presents  a  general  survey  of 
the  basic  mechanisms  for  understanding  local  fallout,  Tftilch  have  been 
studied  In  great  detail.  Since  the  publication  of  this  report,  work  haa 
been  In  progress  toward  establishing  a  better  definition  of  the  problem  pre¬ 
sented  by  these  basic  mechanisms  and  Improving  our  understanding  of  the 
details  Involved.  The  purpose  of  this  paper  Is  to  present  the  results  of 
studies  made  since  the  publication  of  Close-In  Fallout  and  to  bring  out  some 
alight  changes  In  the  concept  of  the  feCLlout  process.  However,  In  order  to 
provide  a  complete  picture,  a  certain  amount  of  the  earlier  background  will 
be  given,  Including  a  brief  resume  of  the  physical  process  of  fEdlout  and 
a  definition  of  the  problem.  Following  this,  each  of  the  iiarameters  which 
enter  Into  the  problem  will  be  tSLken  up  In  detail,  and  the  results  of  studies 
made  over  the  past  18  months  will  be  presented.  The  presentation  will  be  of 
a  conceptuEil  nature  rather  thsui  a  mathematical  formulation. 

Because  of  the  great  amount  of  technical  detail  to  be  presented  In 
this  report,  the  scc^e  of  the  subject  matter  will  be  more  limited  them  In 
Ref.  1.  The  basic  parameters  for  a  8urfeu;e-burst  bomb  will  be  discussed 
thoroughly,  but  the  effect  of  height  above  the  surface  and  of  changing 
surface  conditions  can  be  only  briefly  touched  v^on.  The  extent  to  which 
local  fallout  can  be  determined  will  be  considered,  and  the  prospects  for 
numerical  forecasting  will  be  presented.  Since  this  Is  to  be  a  report  on 
research  done  at  RAND,  it  cannot  be  considered  as  a  general  survey  or  status 
report  on  the  entire  problem  of  fallout .  The  work  of  the  many  other  organs 
izatlons  engaged  in  fallout  reseEurch  has  been  used,  but  no  coo^irebenslve 
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BunsBELry  has  been  inade.  This  report,  then,  is  a  summary  of  the  physical 
processes  and  mathematical  models  of  local  fallout  from  surface— burst  atomic 
bombs.  There  is  no  biological,  and  only  a  minimum  of  radiological,  informa;- 
tlon  contained  herein. 

The  detonation  of^' an  atomic  device  releases  large  amounts  of  energy 
either  through  the  fission  of  heavy  atoms  or  the  fusion  of  light  ones.  'The 
release  of  energy  may  be  accongjllshed  through  fission  only  or  through  a 
combination  of  fission  and  fusion.  The  fiagiucuto  of  the  fission  of  the 
heavy  atcms  are  generally  radioactive  particles.  Another  source  of  radio¬ 
active  particles  or  atoms  Is  the  activity  induced  in  seme  types  of  stable 
atoms  by  the  capture  of  neutrons.  The  nature  of  the  radioactive  products 
therexore  varies  greacj-y,  uepenuAn^  on  whe?  mxjcture  of  fiosion  and  fusion 
in  the  bomb  Eoad  the  elements  which  are  available  nearby  to  capture  neutrons. 
Because  of  the  variation  which  is  possible  in  the  nature  of  the  radloeictlve 
products,  this  report  will  not  deal  with  any  purely  radiological  parameters; 
instead,  the  fraction  of  the  device  will  be  the  primary  Independent  varia¬ 
ble.  The  fractions  which  ai'e  used  and  which  will  be  discussed  further  in 
this  report  are  based  on  measurements  of  aolybdenum^^  (Mo^),  a  radioactive 
element  formed  in  a  predictable  way  by  fission,  which  general  l y  falls  out 
without  fractionation  and  represents  a  very  good  measure  of  the  radioactive 
fraction  of  the  device.  If  the  details  of  the  radiological  con^ionents  of  the 
bomb  and  the  soil  are  known,  it  is  elvayo  possible  to  synthesize  a  proper 
radiological. factor  to  be  multiplied  by  the  fraction  of  the  device,  hence 
generating  either  dose  rate,  integrated  dooeige,  or  any  other  measure  of 


Fractionation  is  the  aeparatloii  of  different  Isotcpes  due  to 
dtffcrencec  in  the  natui'c  of  the  radioactive  precursoro. 
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the  radioactivity. 

The  radioactive  isotopes  which  are  formed  are,  for  the  inost  part,  con¬ 
tained  in  an  extremely  hot  mass  of  air  and  residual  honib  debris*  At  the 
time  of  detonation  most  of  these  Isotopes  will  exist  in  the  form  of  a  cloud 
of  ions.  As  the  cloud  rises  and  cools  the  Isoti^es  form  into  particles  by 
attaching  themselves  to  solid  particles  which  are  swept  through  the  cloud> 
or  by  condensation  into  minute  drops,  or  by  a  combination  of  both  processes* 
Naturally,  the  most  refractory  materials  will  form  into  particles  first  and 
some  of  the  less  refractory  materials  later.  Of  course,  the  noble  gues 
will  not  condense  at  all.  We  have  no  measurements  of  particle  formation 
inside  the  fireball^  thus,  our  knowledge  of  the  actual  method  of  formation 
of  the  particles  is  not  complete.  However,  an  examinatlGn  of  the  resulting 
particles  enables  us  to  estimate  the  distribution  of  their  mass  and  size. 

It  has  also  been  detexmined  that  some  radioactive  elements  which  have  gas¬ 
eous  precursors  do  not  form  particles  early  enough  in  the  process  to  be 
precipitated  in  the  same  way  as  some  of  the  more  refractory  materials. 

Thus,  the  problem  of  the  re fractionation  —  the  atmospheric  separation  of 

some  of  the  radioactive  elements  from  the  main  body  of  the  fallout  —  must 

,  * 
be  borne  in  mind. 

The  particles  which  are  thus  formed  and  the  activity  which  they  carry 
r^resent  one  of  the  necessary  inputs  to  any  fallout  model.  The  informa¬ 
tion  which  we  have  about  the  size,  mass,  and  fall  rates  of  these  radioactive 
particles  and  the  distribution  of  the  radioactivity  with  particle  size  has 
actually  all  been  deduced  from  observations  of  the  fallout  pattern  (the  two 


* 

Strontium  90  and  cesium  137»  two  materials  strongly  fractionating 
seem  to  be  deposited  in  a  way  that  is  different  from  that  of  Mo99, 
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different  ways  lr.i  which  these  deductions  have  been  made  will  be  discussed 
in  detail  later).  This  distribution  represents  one  of  the  la^iortant  Initial 
parameters,  for  it  is  necessary  to  know  what  fraction  of  the  activity  re¬ 
sides  in  particles  of  different  sizes. 

The  tremendous  energy  released  by  a  nuclear  detonation  over  a  small 
space  and  in  a  short  time  causes  a  large  bubble  of  hot  air.  This  hot,  low- 
density  air  has  a  tendency  to  rise  through  the  atmosphere.  As  it  rises  it 
does  work  against  the  atmosphere  and  in  so  doing  is  cooled.  It  has  been 
observed  that  the  outward  gradient  of  ten5>erature,  in  the  hot  bubble  of  air, 
causes  the  central  portion  to  rise  relative  to  the  outside  portion  of  the 
cloud.  This  caxises  an  overturning  of  the  cloud,  and  a  toroidal  or  smoke 
ring  circulation  develops.  Most  of  the  particles  of  radioactive  debris 
end  v?)  in  the  core  of  the  smoke  ring.  This  is  an  extremely  stable  circular- 
tion.  Most  of  the  radioactivity  is  carried  aloft  with  the  cloud,  and  only 
a  emeil  fraction  escapes  from  the  toroidal  circulation  to  be  left  behind 
to  form  a  wake  or  stem.  The  material  thus  carried  toward  forms  the  radio¬ 
active  cloud. 

When  the  energy  which  caused  the  cloud  to  rise  has  been  dissipated  by 
rising  through  and  mixing  with  the  atmosphere,  the  cloud  arrives  at  a 
stabilization  point,  the  height  of  which  depends  very  markedly  on  the  yield 
and,  to  some  extent,  on  the  conditions  in  the  atmosphere.  For  the  larger 
yields,  the  toroidal  circulation  continues  for  a  while  after  stabilization, 
and  this  tends  to  ejqpand  the  radius  of  the  toroid;  thus,  as  the  cloud 
approaches  Its  stabilization  point,  a  continued  horizontal  expansion  is 
noticed.  At  about  five  or  six  minutes  after  detonation,  when  the  cloud  has 
stabilized,  the  initial  space  distribution  of  the  rEidioactivity  is 
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eatabllohed.  It  is  in  the  form  of  a  ring-shaped  body  of  debris  of  particles 
of  variour,  sizes.  In  addition  to  this  ring-sh£Q)ed  body  there  vlU  be  sotae 
trailing  materiel  dovn  along  the  wake  or  stem  of  the  cloud.  The  experiments 
during  Operation  Redwing,  about  which  more  will  be  said  later,  definitely 
show  that  there  is  a  very  small  fraction  of  the  debris  contained  in  the 
stem  of  the  cloud.  The  space  distribution  of  the  radioactivity  represents 
another  lngportant  parameter  of  the  radioactive  fallout.  Those  same  Redwing 
experiments  demonstrate  quite  conclusively  that  the  majority  of  the  radio¬ 
active  debris  is  in  the  lover  portion  of  the  visible  cloud  and  that  it  is 
contained  in  a  ring  arotmd  a  central  portion  which  has  relatively  little  ' 


debris.  .  ■  V 

After  stabilization,  the  tcrcidal  circulation  rapidly  decays,  and  the 
psurtlcles  of  debris  fall  from  that  elevation  to  the  ground.  The  particles 
will  attain  scxne  sort  of  terminal  velocity  and  will,  of  covtrse,  be  trana-’,  . 

■  '  V* 

ported  by  the  wind.  If  it  is  assumed  that  the  vertical  motion  of  the  air 
is  negll2^bl•>  in  con5)arlson  with  the  fall  rate  of  the  particles,  then  the 
vertical  ^-elocities  of  the  particles  can  be  confuted.  Thus,  the  fall  rate 
of  the  particles,  which  will  be  a  function  of  the  size,  shape,  height,  and 

V  ■3. 

density  of  the  particles,  is  an  irportant  and  necessary  conponent  of  the 


rao-j-out  analysis.  ^  • 

It  can  be  shown  that  for  the  particles  which  comprise  the  majority  of  . 

the  fallout>  the  wind  can  be  considered  to  transport  them  with  the  ^eed 

*  »  • 

of  the  wind.  That  is  to  say,  there  is  no  appreciable  lag  in  the  pEirtlclea. 
following  the  wind  flow.  Suffice  it  tc  itate  here  that  the  wind  will  carry 


the  particles  in  the  horizontal  direction,  and  it  is  Imperative  to  know  where 
the  wind  is  going  to  carry  them.  .  ,  ‘ 
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There  are  several  variations  la- the  vd.nd  which  must  be  considered. 

Firet  of  allj  there  are  large  shears  of  the  horizontal  wind  in  the  vertical, 
and  these  should  be  accounted  for.  Furthermore,  some  of  the  particles  take 
quite  a  bit  of  time  to  reach  the  ground,  and  during  this  time  the  wind  patterns 
will  change,  so  the  time  vcLrlation  of  the  winds  must  be  taken  into  account. 
Also,  some  of  the  particles  eure  carried  over  great  distances,  so  we  need 

to  take  into  account  the  variation  of  the  wind  in  the  horizontal.  Thus,  to 

1 

be  precise,  the  ^ind  should  be  considered  as  a  function  of  the  space  dimen¬ 
sions  and  time. 

To  sTmmarlze,  then,  in  order  to  con^cute  or  predict  fallout  it  is 
necessary  to  .know  the  distribution  of  activity  with  particle  size,  the  dis¬ 
tribution  of  the  activity  In  space  at  the  time  of  stabili-AAtlon;  the  fall 
velocity  of  the  particles,  and  the  vind  structure.  In  the  following  sec¬ 
tions  these  parameters  will  be  considered  in  turn.  They  will  be  considered 
chiefly  in  connection  with  surface-burst  bomba  in  the  megaton  reuige. 

Finally,  we  shall  attenqit  a  synthesis  of  the  individual  parameters  into 
the  final  fallout  pattern.  Note  that  there  is  no  necessity  in  this  work 
for  discussing  the  radiological  nature  of  the  debris.  We  speak  of  the 

fraction  of  the  device,  and  we  confute  this  fraction  on  the  basis  of  Mo^^  - 

one  oi  the  oesb  miown  laajLout  con^ionenxs.  Thus,  the  cleanliness  of  the 
weapon,  the  f Issiorr-fusion  ratio,  the  nature  of  fission  products,  or  Induced 
activities,  need  not  concern  mb  In  defining  a  fallout  pattern.  We  will  de¬ 
fine  the  fraction  of  the  device,  and  if  there  Is  any  fxirther  information  as 
to  the  nfiture  of  the  radioactive  material,  this  can  be  incorporated  with 
the  fraction— of— device  numbers  and  el  1  the  customary  radiological  values 
deduced.  Obviously,  since  we  tove  not  considered  radiological  porsmeters. 

We  -vlll  not  concern  oui-nelves  wltli  the  biological  effects  of  radiation. 
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There  have  been  quite  a  few  atten^ts  to  deduce  the  distribution  of 
radioactivity  with  particle  size.  One  of  the  most  usual  methods  of  making 
this  deduction  is  to  plot  the  position  of  various  sized  particles  from 
varioTis  elevations  In  the  clouds  according  to  their  fall  velocity  and  the 
winds,  then  to  measure  the  ground  activity  at  that  point  and  to  assign  the 
fraction  of  activity  as  a  function  of  particle  size  In  accordance  with  the 
amount  on  the  ground  associated  with  that  particle  size.  This  was  essentially 
the  method  used  to  determine  the  particle  size  emd  activity  distribution 
for  the  Jangle  surface  shot  as  presented  in  RAND  report  R-265“AEC.^^^  It 
is  also  the  essence  of  the  Weather  Bureau  procedure  In  determining  the 
activlty^ize  distribution  for  tower  shots. 

Another  more  direct  method  Is  to  simply  measixre  the  activity  on  par¬ 
ticles  that  have  been  sized  after  they  have  been  collected.  This  is  the 
method  used  by  Krey^^^  in  emalyzing  the  activity  and  particle  size  distri¬ 
bution  from  the  collections  on  YAQ-hO  at  the  Zuni  Redwing  shot. 

Figure  1  presents  several  sets  of  results  that  have  been  obtained  In 
the  past.  The  scales  may  be  awkward,  but  they  are  Ideal  for  containing  all 
of  the  data.  The  ordinate  Is  a  logarithmic  scale  and  the  abscissa  Is  the 
Integral  of  the  Gaussian  error  function.  Thus,  a  straight  line  on  this 
particular  graph  would  represent  a  "normal  distribution"  in  the  logarithm 
of  the  particle  size.  The  graph  is  cumulative  so  that  the  fraction  of 
material  on  particles  less  than  the  size  shown  on  the  ordinate  Is  read  on 
the  abscissa.  The  heavy  solid  line  is  the  log  normal  curve  which  is  fitted 
to  the  distribution  taken  from  Ref,  2,  The  circles  are  points  from  the 
distribution  in  Ref .  2,  It  can  be  seen  that  most  of  the  data,  at  least  near 
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the  center  of  the  graph.  Is  clustered  very  closely  around  this  line.  It 
may. perhaps :  be  worthwhile 'to  note  some  of  the, deviations  from  this.  To  the 
upper  lef t  of  the  .strail^t shown  the' ®DL  cxa've  derived  from  Jangle 
surface  and  Linder groLind'  shots  by  analyzing  particle  size. This  shows  a 
much  larger  fraction  of  large  particles.  It  is  believed  that  the  fact  that 
data  from  the  underground  shot  were  used  in  deriving  this  ciorve  has  caused 
it  to  indicate  too  many  large  particles  for  a  surface  burst.  The  short  bit 
of  cujrve  on  the  lower  right  of  the  straight  line  represents  the  average  of 
the  Weather  Bureau  exataiiiatlon  of  the  tower  shots  in  Nevada.  It  is  to 
be  expected  that  the  tower  shots  would  have  a  larger  fraction  of  small 
particles  because  there  is  not  as  much  dirt  taken  into  the  fireball. 
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ing.  The  lower  end  of  the  curve,  that  is,  the  part  dealing  with  the  smaller 
sizes,  shows  an  almost  linear  trend  on  the  log-normal  paper.  However,  for 
particles  larger  than  7**^  microns  the  curve  falls  off,  showing  100  per  cent 
of  the  materled  on  particles  less  than  210  microns  in  radius.  This  is 
probably  caxised  by  the  fact  that  the  larger  particles  fall  to  earth  long 
before  they  coLild  reach  the  collection  station.  Thus,  the  curve  I’epresent- 
Ing  Krey's  dAta  is  based  on  a  biased  sanqple.  The  same  may  be  true  for  the 
Ref.  2  curve  wnicn  lea-Ls  uxa  ah  a  similar  manner  but  wnich  reacnes  nigner 
values. 


A  log-normal  curve  appears  to  be  a  good  fit  to  the  hLilk  of  the  data. 
The  Weather  Bureau  curve  for  tower  shots,  which  appears  at  the  lower  right 
of  the  SLxrf ace— burst  data,  can  be  discounted  because  they  are  not  surface 
bursts.  The  NRDL  cui-ye,  at  the  upper  left  of  tJie  main  body  of  data,  can 
be  discoLinted  becuLise  of  the  inclusion  of  data  from  the  \inderground  shot. 
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The  efficiency  of  large  particles  in  Krey'a  curve  and  the  estimate  presented 
in  Ref.  2  can  he  attributed  to  the  removal  of  the  large  particles  in  such 
short  time  that  they  do  not  enter  into  the  san^illng  procedures. 

The  mantier  in  which  the  curves  presented  in  Fig.  1  have  been  obtained 
is  more  or  less  biased.  There  has  been  no  method  found  to  date  for  getting 
a  completely  vinbiased  earple  of  the  particle  size  from  the  original  cloud 
and  the  prospect  of  obtaining  unbiased  measurements  is  not  bright^  because 
of  the  tremendous  problems  Involved  in  sampling,  sizing,  and  measuring  the 
radioactive  particles. 

The  work  done  to  date  has  servea  to  indicate  that  a  log-normal  distri¬ 
bution  appears  to  be  a  reasonable  choice  for  the  distribution  of  activity 
with  particle  size.  There  has  been  some  evidence  that  there  is  a  different 
distribution  of  activity  with  particle  size  at  different  elevations;  in 
other  words,  the  distribution  of  activity  with  particle  size  is  not  inde¬ 
pendent  of  height.  However,  for  the  purposes  of  this  study,  the  assumption 
was  made  that  the  activity  distribution  with  particle  size  was  completely 
Indcpendeht  of  the  activity  distribution  Tfith  height.  A  review  of  aome  of 
ths  attempts  to  separate  the  size  distributions  according  to  height  indicates 
that  the  assumption  of  independence  will  not  cause  any  serious  difficulty. 

As  a  matter  of  fact,  none  of  the  observations  are  sufficiently  accurate  to 
warrant  the  conclusion  of  differing  particle  size-activity  distributions 
with  height. 

Since  the  distribution  appears  to  be  a  log-normal,  and  since  the  data 
studied  BO  far  has  considerable  bias,  it  was  decided  to  attempt  to  recon¬ 
struct  a  log-normal  particle  size-activity  distribution  which  would  best 
fit  the  Zunl  reaulto  on  Operation  Redwing.  The  basic  assumptions  of  this 
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approach  are  that  the  dletrlhutloa  is  log-aomal,  that  it  Is  Indepeadeat  of 
elevation  (as  mentioned  above),  that  the  distribution  of  activity  with  heighjb 
is  known,  and  that  the  winds  are  known  as  functions  of  space  and  time.  The 
excellent  series  of  weather  maps  which  were  analyzed  by  the  Joint  Task  Force  ^ 
Weather  Group  for  the  Redwing  Tests^^^  was  used  to  construct  time-varying  and 
space-varying  wind  patterns.  The  distribution  of  the  activity  with  height 
was  taken  from  the  rocket  resTilts  as  reported  in  "Fallout  Studies  During 
Operation  Redwing.  Using  the  method  of  computation  as  outlined 
in  Ref.  8,  five  different  log-normal  distributions*  were  chosen  and  each 
of  them  was  used  to  make  a  computation  of  the  activity  at  each  of  five 
different  measuring  stations*  For  each  station  a  plot  was  made  on  axes 
which  represented  the  mean  of  the  logarithm  of  particle  size  for  the  ordinate 
and  the  standard  deviation  of  the  logarithm  of  particle  size  for  the  abscissa* 
Each  of  the  five  distributions  on  this  graph  are  therefore  r^resented  by  a 
point.  At  each  of  these  five  points  the  confuted  fraction  of  the  device  tot 
that  particular  meeisurlng  station  was  entered.  The  five  values  of  fractions 
Were  then  used  to  STibjectlvely  construct  lines  of  equal  fraction  of  device 
as  a  function  of  mean  of  the  logarithm  of  particle  size  and  standard  deviar~ 
tion  of  the  logarithm  of  particle  size.  The  measured  value  (flee  Ref.  9)  vaS 
then  used  to  find  the  line  of  possible  means  and  standard  deviations  for 


The  log-normal  function  Is  simply  the  GaussiEm  distribution  in  the 
logarithm  of  the  variable.  If  x  -  log  r  and  the  distribution  of  x  is  normal, 
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This  is  a  two-parameter  distribution  in  which  the  distribution  is  conroletely 
determined  by  the  means  of  the  logarithm  of  the  particle  radius  eind  the 
standard  deviation  of  the  logarithm  of  the  particle  radius. 
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each  of  the  stations.  On  Fig.  2  the  observed  range  of  each  of  the  five 
stations  Is  siQ)eriii?)osed.  It  may  be  noted  that  there  are  no  lines  for  the 
results  from  YAC>-59»  This  station  was  at  an  extreme  edge  of  the  fallout 
pattern  and  the  fallout  did  not  arrive  here  until  almost  2k  hours  after  the 
event.  With  this  long  delay  time  the  wind  errors  accumulate  to  such  an 
extent  that  no  confidence  can  be  placed  In  the  wind  plot  (see  Sec.  V). 
Therefore,  lack  of  fit  of  this  one  station  does  not  seriously  discredit 
the  results. 

The  other  foLjr  stations  show  that  the  mean  and  standard  deviation 
which  best  fit  all  of  the  stations  must  be  within  the  area  that  has  been 
chosen.  Subjective  best  estimates  for  the  mean  of  the  logarithm  of  the 
radius  and  the  standard  deviation  of  the  logarithm  of  the  radius  are  I.65 
and  0.30  respectively.  This  particle  size  distribution  was  then  used  in 
the  model  with  the  wind  from  the  Tewa  shot  from  Operation  Redwing.  The 
results  that  this  particle  size-activity  distribution  gave  on  the  independr- 
ent  test  with  the  Tewa  shot  were  most  gratifying,  A  plot  of  the  con^iuted 
versus  the  measured  fraction  is  shown  in  Fig.  3» 

Because  of  the  known  uncertainty  in  the  particle  size  distribution,  it 
is  necessary  to  know  something  about  the  effect  of  any  changes  in  the  dis¬ 
tribution  on  the  fallout  pattern.  It  is  apparent  that  a  large  fraction  on 
small  particles  will  produce  an  extensive  pattern.  A  large  fraction  on 
large  particles  will  produce  a  less  extensive  pattern  but  one  of  greater 
intensity.  For  a  log-nonnal  distribution  with  a  given  standard  deviation, 
an  increase  in  the  mean  will  place  more  activity  on  leirger  particles  which 
sure  deposited  near  ground  zero  and  less  on  particles  which  are,  for  the 
most  part,  too  small  to  appear  in  the  local  fallout  in  s^  event..  Thus 
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Fig.  2  —  Subjectively  determined  combinations  of  log  r  and 
^log  I-  'which  would  fit  within  the  stondard  error  of 
observation  of  fraction  down  at  4  observing 

stations 
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Observed  fraction  x  10' 

Calculated  versus  observed  fraction  of  device  per  square  kft  for  the  Tewd  shot— Operation  Redwing 
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the  changes  In  the  vaGLus  of  the  mean  of  the  logeorlthm  produce  large  changes 
in  the  pattern  near  ground  zero.  For  example,  a  change  in  the  mean  of  the 
log  of  r  from  1.47  to  I.67  increaued  the  fraction  at  the  statlona  within 
20  mi  of  ground  zero  by  factors  of  from  4  to  10.  At  more  distant  statlona, 
where  particles  from  the  center  of  the  distribution  cause  the  fallout,  the 
shift  of  the  mean  will  cause  leas  change  in  the  fraction  in  any  particle  size 
range,  therefore  a  shift  in  the  mean  of  the  magnitude  quoted  above  produces 
an  Increase  of  a  factor  of  only  1.5  for  a  station  — ^  approximately  60  mi 
from  ground  zero.  On  the  distant  edge  of  the  local  fallout  pattern  only  a 
fine  rain  of  small  particles  occurs  and  the  fraction  is  quite  insensitive 
to  change  in  the  mean. 

For  a  given  mean  value,  an  increase  in  the  standard  deviation  of  the 

logHiozmial  will  cause  an  increase  in  the  fraction  of  particles  at  both  ends 

of  the  distribution.  Since  the  small  end  of  the  particle  spectrum  does  not 

contribute  much  to  the  fallout,  the  change  on  the  smaJl  end  is  not  noticed. 

The  depletion  of  the  center  of  the  distribution  shows  i^^ae  a  decrease  in 

fraction  at  intermediate  distances  with  increasing  a j  the  augmentation  of 

the  large  end  shows  up  as  an  increase  in  fractions  at  very  close  stations 

with  increasing  a.  Thus  an  increase  in  o,  from  .27  to  .57  decreases 

log  r 

the  fraction  slightly  in  the  area  where  the  center  of  the  distribution  falls 
and  doubles  or  triples  the  fraction  in  the  region  close  to  ground  .zero. 

We  believe  that  the  distribution  of  activity  with  particle  size  is 

determined  more  by  the  conditions  of  the  burst  - the  tenperature  obtained 

and  the  nearness  to  the  surface  — than  by  the  nature  of  the  underlying 
surface.  Until  further  evidence  to  the  contrary  can  be  presented,  we  will 
assume  that  the  distribution  of  activity  with  p^ticle  size  for  surface- 


CONFIDENTIAL 


rm-255*^ 

2-25-59 

16 


CONFIDENTIAL 

FORMERLY  RESTRICTED  DATA 

ATOMIC  EHCRCr  ACT— 1954 


burst  veapons  will  b«  that  of  a  Ibg-normal  vlth  a  mean  of  1*65  and 
standard  deviation  of  O.JO, 
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The  second  parameter  vhlch  must  be  taken  iqi  in  some  detail  is  the  die- 
tributlon  of  the  radioactive  peirtlcles  in  space  at  the  time  of  atabilizet- 
tlon  of  the  cloud.  The  study  of  the  space  distribution  of  activity  has  a 
long  history  and,  in  the  early  days,  aome  very  conflicting  reports.  The 
early  reports,  that  la,  the  reports  of  all  of  the  tests  t?)  to  the  Redwing 
series,  were  based  on  observations  of  the  visual  cloud.  For  some  of  the 
smaller  yields,  the  toroidal  circulation  was  not  apparent,  and  it  was  there¬ 
fore  thought  that  there  was  a  Jet  of  material  blasted  into  the  air,  and 
that  there  was  considerable  mater iad  at  all  levels.  However,  after  some  of 
the  larger  shots  in  the  Pacific,  when  the  toroidal  circulation  became 
evident  and  after  some  of  the  reconstructions  which  were  done  at  the  Weather 
Bureau^^^  and  at  RAND, it  was  realized  that  the  activity  was  largely 
confined  to  the  mushroom,  Euid.  that  there  should  be  a  concentration  in  the 
core  of  the  torus.  There  was,  however,  no  evidence  as  to  the  ratio  of 
activity  in  the  mushroom  or  top  part  of  the  cloud  to  that  which  was  in  the 
stem. 

Based  on  a  reconstruction  from  the  fallout  pattern.  Ref.  2  indicated 
90  per  cent  in  the  mushroom  and  10  per  cent  in  the  stem.  The  dimensions 
of  the  cloud,  which  were  used  in  all  of  the  early  reports,  are  those 
observed,  cilthough  it  was  recognized  that  water  vapor  played  a  great  role 
in  defining  the  limits  of  the  visible  clotid,  particularly  near  the  stabi¬ 
lization  time.  Kellogg^^^^  has  reported  in  great  detail  on  the  early 
studies  of  dimensions  of  the  cloud  as  well  as  the  means  for  measuring  and 
correcting  for  such  things  as  atmospheric  winds  and  opticeil  problems . 

One  of  the  problems  which  arose  with  the  Castle-Bravo  test  was  the 
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tremendous  horizontal  expansion  of  the  cloud  Just  at  stabilization.  It  was 
pointed  out  in  Ref.  11  that  if  the  radioactivity  was  indeed  spresid  as  far 
as  the  visible  cloud  and  if  constant  concentration  with  radius  was  assumedi 
it  would  be  difficult  to  explain  the  high  concentration  of  fallout  that  was 
observed  from  the  Castle-Bravo  shot.  Therefore,  the  idea  of  an  effective 
radiological  diameter  which  was  somewhat  smaller  than  the  final  visible 
diameter  was  Introduced. 

The  observations  of  the  toroid  and  studies  of  smoke  rings  suggested 

that  most  of  the  debris  should  be  concentrated  in  a  doughnut-shaped  space 

in  the  clovid.  It  was  not  until  rockets  measuring  radioactivity  were  shot 

through  the  clouds  during  severed  of  the  Redwing  tests,  however,  that  this 

(7) 

hypothesis  was  adequately  confirmed. ' ' '  These  shots  also  indicated  that 

the  amount  of  radioactivity  in  the  outer  reaches  of  the  visible  cloud  was 

Indeed  extremely  small,  and  that  the  effective  radiological  diameter  was 

much  smaller  than  that  of  the  visible  cloud. 

The  radiological  measurements  which  were  made  from  Redwing  rockets 

sinply  reported  the  flux  of  gamma  rays  or  the  roentgen  dose  encountered  by 

rockets  as  they  flew  along  trajectories  through  the  clouds.  In  "Fsdlout 

(7) 

Studies  During  Operation  Redwing,"'’'  the  first  reports  of  the  results  of 
these  rocket  measurements  were  given.  The  only  set  of  measurements  which 
were  of  particular  interest  in  this  project  were  those  of  the  Zuni  cloud. 
Other  clouds  were  measured  but  they  were  formed  by  alrburst  or  barge  shots; 
also,  in  some  cases,  too  few  rockets  penetrated  the  clpvid  to  make  the 
analysis  usable. 

The  basic  rocket  data  was  greatly  smoothed  and  idealized  to  obtain  the 
desired  distribution  from  the  Zuni  data.  These  snwothed  YfOLues^were.  numerrr 
Ically  integrated  over  the  volume  swept  out  by  the  figures  of  revolution 
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formed  by  rotating  the  contour  lines  about  the  axis.  Each  grid  valvie  vas 
then  divided  by  the  Integrated  value  to  yield  the  fraction  of  the  device 
per  cubic  kft  at  each  of  the  grid  Intersections.  These  values  Vere  then 
Integrated  at  each  height  over  the  area  of  the  circle  defining  the  cloud  at 
that  level  In  order  to  arrive  at  the  distribution  of  activity  as  a  function 
of  height.  Flgwe  4  shows  this  distribution  In  fraction  per  kft.  The  dls-* 
trlbution  shown  In  Fig.  4  is  a  very  much  idealized  representation  of  the 
fraction  of  device  per  kft  for  the  Ztinl  shot  of  Operation  Redwing. 

In  order  to  determine  the  radial  distribution,  the  values  of  concerr* 

tratlon  were  integrated  over  height  for  various  values  of  radius.  The 

2 

resulting  distribution  of  activity  per  kft  as  a  function  of  radius  is 
shown  in  Fig.  5>  Again,  this  Is  an  Idealized  representation  of  the  Zunl 
shot,  but  it  does  show  the  reality  of  the  concentration  in  the  toroid. 

The  pee*  concentration  close  to  the  center  and  the  small  values  at  great 
distances  demonstrates  why  the  combination  of  the  visible  cloud  reidlus  with 
a  constant  radial  distribution  failed  to  produce  sufficiently  Intense 
radiation  on  the  Castle-Bravo  shot. 

Before  going  on  with  the  problem  of  scaling  these  results  to  different 
yields,  it  is  necessary  to  discuss  the  sensitivity  of  the  fallout  pattern 
to  variations  of  these  particular  distributions.  The  radial  distribution 
will  be  considered  first  because  it  has  little  effect  on  the  pattern. 
Consider  a  single  layer  of  particles  of  varying  sizes.  As  the  particles 
fall  to  the  ground  at  different  rates,  the  overlapping  of  the  rings  will 
tend  to  obscure  the  low  concentrations  in  the  middle  of  the  initial  distri¬ 
bution.  In  order  to  demonstrate  this  effect  in  on  actual  pattern,  the 
computing  model  was  used  to  coapute,  with  a  very  simple  wind  structure,  a 
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Fig.  5  —  Smoothed  distribution  of  activity  with 
radial  distance  from  cloud  center  from 
Zuni  rocket  data 
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pair  of  patterns^  one  based  on  the  distribution  shown  in  Pig.  5  *nd  the 
other  on  a  constant  concentration  to  an  effective  radiological  reuilus. 

Figure  6  shows  a  conparlson  of  the  cross  section  of  the  resultsmt 
fedlout  pattern  for  two  radial  distributions  at  a  point  downwind  which  is 
representative  of  all  downwind  cross-sections.  It  can  be  seen  that  the 
correct  radial  distribution  yields  a  more  peaked  cross-eectlon  with  a 
slightly  wider  base,  Indicating  a  spread  of  small  fractions  over  a  wider 
area.  As  will  be  seen  later,  the  edges  of  the  pattern  cannot  be  precisely 
determined  in  any  event,  so  that  the  use  of  an  effective  radiologicsLl  radlUs 
with  a  constant  concentration  is  a  good  method  of  approximating  the  radial 
distribution. 

IThe  effect  of  cnauges  in  the  elevation  of  the  activity  peak  shown  in 
Pig.  4  was  also  checked  by  machine  computations.  It  would  have  been  deslrm- 
ble  to  check  the  entire  course  of  the  distribution  but  there  were  not  enough 
ground  measurements  to  provide  a  reasonable  test.  In  order  to  test  the 
effect  of  changes  in  height  of  the  distribution,  the  shape  of  the  distrlbvr' 
tlon  cvirve  was  preserved.  Calcvslations  were  made  for  the  fraction  down  for 
the  Zuni  shot  from  the  distribution  as  shown  In  Fig.  4,  Emd  also  from  a 
distribution  with  the  height  of  the  peak  concentration  raised  by  5000  ft. 

The  variation  of  fallout  at  the  close-in  stations  was  practically  unchsinged, 
but  at  YAQ-40,  further  out,  there  was  a  marked  increase.  This  is  due  to 
the  fswt  that  material  from  the  peak  of  the  activity-height  distribution 
iras  made  to  reach  yAQ-4o  when  the  peak  height  was  raised.  Since  the  original 
height  distribution  fit  YAO-4o  very  well  over  a  rather  large  range  of 
actlvity^slze  distribution,  it  may  be  concluded  that  the  activity  distribution, 
derived  from  the  rockets  as  shewn  in  Fig.  4,  was  about  correct. 
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Fig.  6  —  Cross  section,  80  miles  downwind  of 
ground  zero,  for  two  radial  distributions 

Washington  D.C.  mean  winter  winds 
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Having  determined  distributions  of  activity  vith  height  from  the  Zunl 
shot^  it  is  necessary  to  provide  some  method  for  scaling  these  distributions 
to  other  yields.  For  such  a  scaling,  recourse  must  be  made  to  visual  obser¬ 
vations  because  the  radiological  measurements  are  insufficient.  The  varia¬ 
tion  cf  the  center  of  a  cloiid  in  the  tropics,  as  deduced  by  Kellogg, 
is  shown  in  Fig.  7  together  with  the  height  of  the  top  euid  bottom  of  the 
Redwing  clouds  as  reported  at  seven  minutes  in  Cloud  Photography.  The 

agreement  between  Kellogg's  analysis  and  the  observations  is  extremely  good. 
Rough  estimates  of  the  height  of  the  maximum  of  activity  for  Zuni,  Cherokee, 
and  Navajo  shots  are  shown  by  the  large  circles  in  Fig.  7.  The  best  scaling 
that  can  be  deduced  from  this  meagre  data  is  that  the  peak  height  of 
activity  is  coincident  with  the  base  of  the  visible  cloud.  The  line  of 
best  fit  to  the  base  of  the  visible  cloud  can  be  expressed  as 

H{kft)  =  42 


for  yields  greater  than  500  KT  in  tropical  atmospheres.  For  other  atmos¬ 
pheric  conditions,  it  is  recommended  that  Kellogg's  curves  be  used  with  a 
correction  to  the  base  of  the  visible  cloud. 

In  order  to  scale  the  radial  distribution,  the  concept  of  effective 
radiological  radius  with  a  constant  concentration  will  be  used.  An  effect 

tive  radiological  radius  of  85  kft  with  a  concentration  of  4.4  x  10~^  parts 
2 

per  kft  was  chosen  as  a  fair  representation  of  the  Zuni  distribution 
shown  in  Fig.  5*  A  plot  of  cloud  radius  at  the  time  of  stabilization 
at  end  of  growth  (see  Fig.  8)  was  taken  from  Ref.  10.  To  this  plot  was 
added  the  seven-minute  visual  radius,  as  taken  from  Ref.  12.  The  single 
point  of  the  Zuni  effective  radiological  radius  was  plotted  eilso.  The 
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scaling  of  the  effective  radiological  radius  vb,s  then  ta.-ken  to  be  a  line 
through  85  kft  at  3. If  MT,  vhich  is  congruent  vith  the  other  curves  of  visi¬ 
ble  radius.  This  may  be  ejqpressed  in  the  range  above  300  KT  as 


R  (kft) 45  (MT) 

6 


These  scaling  lavs  should  be  good  approximations  for  most  cases.  It  must  be 
eiqphaslzed,  however,  that  the  equations  are  derived  for  tropical  atmospheres. 
By  a  JudlcioTis  subjective  correction  for  atmospheric  conditions  following 
the  guide  lines  of  the  Machtar-Kellogg  theory  of  the  rise  of  the  cloud,  reet- 
sonable  resolutions  of  the  radiological  scaling  may  be  made. 
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At  the  time  of  stabilization,  the  activity  is  distributed  over  a  region 
of  space  and  most  of  it  is  carried  by  small,  solid  particles.  A  study  of 
these  particles  indicates  that  they  are  irregular  shapes  with  a  density  of 
about  2.5.  The  density  is  apparently  partly  a  function  of  the  material  of 
vhlcb  they  are  con^iosed  and  the  nature  of  their  formation.  The  value  of 
2.5  has  been  chosen  as  being  very  close  to  the  avereige  value  of  particles 
formed  from  both  Nevada  sand  and  Pacific  cored.. 

As  the  circulation  of  the  cloud  dies  out,  gravity  will  impart  a  dovna- 
ward  force  which  will  rapidly  be  balanced  by  the  viscous  force  of  the  air. 

The  resultant  terminal  velocity  of  the  particles  will  cause  them  to  settle 
to  the  earth  at  a  rate  dependent  i?)on  their  size.  It  is  obvious  that  ■rariar- 
tions  in  density  and  shape  will  cavise  variations  in  the  fall  rate,  but,  for 
a  first  approximation,  it  can  be  assumed  that  the  particles  are  spheres  of 
density  2.5. 

There  are  three  different  regimes  which  govern  the  terminal  velocity 
of  the  particles  in  still  air.  The  extremely  small  particles  in  the  high 
atmosphere  may  be  the  same  order  of  magnitude  as  the  mean  free  paths  of  the 
molecLiles  of  airj  therefore  the  corrections  for  mean  free  paths  must  be  used. 
The  particles  in  the  range  from  about  15fA  radius  to  about  25ti  radius  are 
large,  conqpared  to  the  mean  free  path,  and  yet  have  sufficiently  small 
velocities  to  permit  the  use  of  Stokes'  approximation.  The  larger  particles 
must  be  treated  by  the  aerodynamic  equations  involving  the  dreig  coefficient, 
which  is  a  function  of  the  Reynolds  nLomber,  which,  in  turn,  is  a  function 
of  the  speed  of  the  particle  and  the  particle  size . 

In  the  Stokes  law  range,  the  familiar  equation  for  fall  velocity  is 
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„  2 

„  ^  2  r 
8  9  n 


where  r  is  the  particle  radius,  p  la  the  particle- density,  g  is  the  acceler¬ 
ation  of  gravity,  and  t]  is  the  coefficient  of  viscosity  of  the  air.  For 
the  oTnan  particles  this  must  be  modified  by  the  effect  of  the  mean  free 

(13) 

paths.  The  correction  factor  for  mean  free  paths,  according  to  Kennard, ' 


M  -  1  + 


25  +  0.41  exp 


0.88  r 


where  M  is  the  factor  which  must  be  multiplied  by  W  to  obtain  the  true  ter— 
thTtiaI  velocity  and  I,  is  the  mean  free  path  of  the  air  molecules.  For  the 
larger  particles  the  viscous  force  on  the  particles  is  given  by 


Fj,  -  I  P^TT  Cp 


Trtiere  Fp  is  the  viscous  force  on  the  particles,  is  the  large  particle 
velocity,  p^  is  the  air  density  and  Cp  is  the  drag  coefficient.  The  relation 
of  the  drag  coefficient  to  the  Reynolds  number  for  smooth  spheres  lb  shown 
in  Fig.  9.  The  Reynolds  number,  in  turn,  is  a  function  of  W^,  r,  p^^,  and  tj 

2  p  r  W 
«  r 


Dy  equatiuii  the  drag  force,  i>-  ,  to  the  gravitational  force,  F  ,  where 

u  s 


F  »  ^  7r  0  e 
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an  equation  Involving  the  speedy  the  radius^  the  drag  coefficient^ 

and  various  constants  is  obtained.  Figure  9  relates  Cj^  to  the  Reynolds 

number  and  the  definition  of  the  Reynolds  number  relates  it  to  the  siae 

and  the  speed.  Because  the  relationship  defined  by  Pig.  9  ts  not  amenable 

to  a  mathematical  analysis^  the  confutations  of  W  vere  performed  by  an 

(ik) 

iterative  scheme  outlined  by  Langmuir.'  '  Figure  10  shows  the  vertical 
velocities  of  a  selection  of  particles  representing  the  full  range  of  ter¬ 
minal  velocity  regimes.  The  New  Standard  Atmosphere was  used  for  the 
atmospheric  parameters.  The  correction  factor  for  the  0.5p  particle  at 
loo  kft,  is  ^10.  At  the  other  end  of  the  graph,  the  aerodynamic  fall  of  a 
60Op  particle  is  a  factor  of  10  smaller  than  rate  of  fall  of  the  Stokes  lav* 

There  is  considerable  uncertainty  about  using  the  drag  coefficients  fojf 
spheres  to  calculate  the  terminal  velocities  of  fallout  particles.  In  order 
to  assess  the  possible  error,  an  experiment  was  performed  to  measure  the 

'ft—*. 

drag  coefficient  of  simulated  fallout  particles  as  a  function  of  the  Reynolds 
number  (Ref,  17).  Calibration  experiments  on  spheres  indicated  that  tech¬ 
nique  was  accurate  to  within  about  five  per  cent.  The  results  on  the  fall¬ 
out  particle  models  showed  only  slight  deviations  from  this  ^here  approxi¬ 
mation  except  for  four  particles  which  were  several  times  larger  in  one 
dimension  than  in  the  other  two.  Since  most  fallvout  particles  do  not  deviate 
greatly  from  a  spherical  form,  the  drag  approximation  could  Introduce  an 
error  of  less  than  10  per  cent,  a  small  error,  compared  with  some  of  the 
other  uncertainties  in  the  data. 

The  terminal  velocities  of  the  particles  Eire  only  part  of  the  total 
vertical  transport.  The  effects  of  the  vertical  motions  of  the  atmosphere 
must  be  Included  to  determine  the  true  transport  relative  to  the  earth.  This 
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effect  is  not  sin^ily  an  additive  one  becaxise  the  relative  motion  of  the  air 
past  the  particles  produces  accelerations  vhich  are  not  necessarily  linear. 

In  order  to  make  a  prellmin£u*y  estimate  of  the  effect  of  vertlceJ.  motldn 
df  the  atmosphere ;  the  Smoluchovskl  equation  may  be  used 


dP 


+  c 


vhere  P  is  the  probability  that  a  particle  originating  at  a  height  will 
arrive  at  a  height  z  at  some  time  t.  K  is  the  coefficient  of  diffusion  and 
<s  is  the  terminal  velocity  of  the  particle.  The  assunptions  >dilch  lead  to 
thiti  equation  are  that  the  air  motions  affecting  the  particles  Ere  conqpletely 
reindom,  euad  that  the  force  of  attraction  is  constant. 

As  can  be  seen,  this  is  merely  eui  extension  of  the  diffusion  equation 
to  account  for  the  gravitational  attraction.  This  eqtxatlon  can  be  solved 
with  the  boundary  conditions;  (1)  that  the  particle  is  at  height  z  ■  z 

o 

With  probability  1  at  t  »  Oj  and  (2)  that  the  particles  stick  to  the  ground 
when  they  reach  z  =  0.  The  equation  for  the  probability  that  a  particle 
Will  reach  the  ground  between  t  Euid  t  +  dt  is  then 

P(0,t)  =  (4Trt  (z^/t)  exp  “ 

This  equation  was  numerically  integrated  for  z^  =  17  km,  K  =  1000  m^/mln  and 
2 

4000  m  /min  for  values  of  C  *  2  m/min  and  12.5  m/mln.  The  resulting  curves 
are  shown  in  Fig.  11.  The  velocity  Value*' correspond,  roughly  to  a  lOu 
particle  and  a  25;i  particle  respectively. 

Figure  11  should  not  be  considered  as  a  quantitative  estimate  of  the 
probability  of  the  particle  reaching  the  ground  in  a  given  time,  but  it  does 


2 
z 

P.  _  __0 
m  2  D 


Cz 


c^t 

^Fd 


CONFIDENTIAL 


4000 


CONFIDENTIAL 

FORMERLY  RESTRICTED  DATA  '  35 

ATOMIC  ENERGY  ACT— (954 


CONFIDENTIAL 


rm-2354 

a-25-59 

36 


CONFIDENTIAL 

FORMERLY  RESTRICTED  DATA 

ATOMIC  ENERGY  ACT— 1954 


provide  a  qualitative  estimate  of  the  variation  In  position  vhlch  must  be 
Expected  due  to  the  vertical  conponent  of  turbulent  fluctuations*  Thu8> 
the  ccn^iutatlon  of  the  position  of  a  lOu  particle  which  neglects  the  verti¬ 
cal  turbulence  Is  almost  certain  to  be  wrong.  On  the  other  hand,  the  prob¬ 
ability  Is  large  that  the  2%  particle  will  land  within  a  few  hours  of  the 
cooQiuted  time.  For  larger  particles,  of  course,  the  spread  will  be  even 
less. 

For  very  large-sceile  vertical  motions  of  the  atmosphere  it  would  be 
entirely  proper  to  add  the  vertical  velocity  to  the  terminal  velocity  In 
order  to  determine  the  vertical  transport.  This  procedure,  however,  would 
be  of  little  practical  use,  because  large-^scale  motions  have  small,  magnitudes 
relative  to  the  terminal  velocity  of  large  particles,  and  wnall  meignltude 
with  respect  to  the  small-scale  uncertainties  for  the  amal 1  particl.es. 

These  vertical  transport  effects  define  a  limit  to  the  utility  of  the  direct 
cos^utatlon  of  fallout.  For  particles  which  take  more  than  about  12  hours 
to  fall  through  still  air,  the  error  caused  by  vertical  turbulence  is  large 
enoTJgh  to  ccm^iletely  Invalidate  the  confuting  process.  It  may  be  possible 
to  develop  procedures  for  statistically  including  turbulent  effects  on  the 
particles  to  arrive  at  expected  or  most  probable  values  of  deposition,  but 
the  crude  approach  presented  here  Is  not  sufficiently  precise  to  attenpt 
such  procedures  at  this  time. 
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V.  HORIZONTAL  TRANSPORT 

As  the  particles  fall  from  their  positions  in  the  stabilized  cloud, 
they  are  transported  horizontally  by  the  wind.  The  force  that  the  horlzon- 
teJ.  coD^onents  of  the  wind  exert  on  the  particles  is  the  same  sort  of  viscous 
or  aerodynamic  force  as  that  discussed  in  the  previous  section.  In  the 
absence  of  a  driving  force,  however,  the  relative  velocity  is  rapidly  re~ 
duced  to  zero  and  the  particles  travel  essentially  with  the  wind.  This  canr- 
not  be  precisely  correct,  because  changes  In  the  horizontal  wind  speed  will 
require  accelerations  of  the  particles.  The  lag  time  of  even  the  largest 
particles,  however.  Is  apparently  small  cpnpared  with  the  rate  at  which  the 
particles  esqierience  velocity  changes;  therefore  it  may  safely  be  assumed 
that  the  horizontal  wind  velocity  is  equal  to  the  horlzontEd  particle 
velocity. 

If  \v(x,y,h,t)  is  the  horizontal  air  velocity  following  the  particle, 
then  the  horizontal  displacement,  D,  is  given  by 

T 

»  -  /  \v(x,y,h,t)  dt 
o 

If  it  is  assumed  that  the  vertical  velocity  is  determined  by  the  terminal 
velocity,  the 

W(r,h)  .  § 

and  the  transport  equation  may  be  written  as 


dh 


®(r,h)  =  /  \V(x,y,h,t) 


hi 


where  it  has  been  indicated  that  D  is  a  function  of  the  particle  size  and 
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the  initial  height  of  the  particle.  The  as8un?)tions  regarding  the  vertical 

velocity  W  have  already  been  discussed. 

The  problem  posed  by  the  nature  of  the  horizontal  vind  must  be  considered. 

In  moat  of  the  early  vork,  the  wind  vas  considered  as  variable  only  In  the 

vertical.  The  effect  of  time  and  space  variations  were  included  In  some  of 

(5) 

the  reconstructions  of  initial  cloud  parameters  by  the  Weather  Bureau.'  ' 

The  model  presented  in  A  New  Model  For  Ffillout  Calculations provides 
for  an  interpolation,  in  time  and  space,  from  observed  winds,  to  estimate 
the  vind  on  the  particle  at  a  rather  close  grid  of  points  in  time  and  spause. 
Any  deviations  of  the  actual  wind  on  the  particle  from  the  Interpolated  wind 
Will  cause  errors  in  the  coirputed  position. 

An  approximate  analysis  of  the  errors  in  displacement  were  made  to 
Indicate  the  magnitude  of  this  error.  The  details  of  the  anELlyslB  are  pre¬ 
sented  in  Ref.  l8.  Table  1  shows  the  ratio  of  the  circular  standard  error 
of  the  displacement  vector,  o^,  to  the  assumed  circular  standard  error  of 
the  deviations  of  the  actual  wind  from  the  Interpolated  wind,  for  a  few 
particles.  This  table  shows,  as  did  Fig.  11,  in  Sec.  IV,  that  the  smaller 
particles  cure  so  subject  to  small-scale  errors  that  precise  location  is 
Inposslble. 

Table  1 

RATIO  Ci;'  Op/o^  AS  A  FUNCTION  OF  DnITIAL  HEIGHT  AND  PARTICLE  RADIUS 


Height 

1  Radius  (u) 

(I^t) 
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.747 
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.063 

20 

4.160 
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30 
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1.196 

.375 

.095 
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5.633 

1.321 

.412 

.105 
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6.578 

1.506 

.464 
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.491 
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effect  of  the  magnitude  of  the  "irind  on  h  fallout  pattern  is  an 
in^iortant  consideration.  In  the  case  of  no  wind,  the  pattern  vould  consist 
of  the  projection  of  the  vortex  ring  i^on  the  ground,  Bynmetrlcailly  around 
ground  zero.  This  vould  be  a  relatively  small  pattern  of  high  concentration. 

On  the  other  hand,  a  strong  vlnd  wovild  carry  small  particles  greater  dls-* 
tances  than  the  large  particles  emd  vould  therefore  produce  a  more  extensive 
pattern,  but  one  of  lesser  Intensity.  In  order  to  eramine  this  effect  the 

(19) 

"single  vind  model,"  developed  in  Simplified  Model  For  Fallout  Calcvilations^ 
and  briefly  outlined  in  Sec.  VI,  vas  employed.  The  results  are  shown  in 
7igs.  12-lli'.  Figure  12  shovs  how  the  value  of  the  maximum  fraction  of  the 
debris  (roughly  the  same  as  the  maximum  "ihot  spot" )  decreases  with  increae-* 
ing  wind  velocity.  Figure  13  shows  the  distance  from  ground  zero  at  ^Ich 
this  maximum  fraction  is  found.  The  dashed  portion  of  the  line  on  Fig.  13 
indicates  that  the  maximum  over  this  remge  of  winds  is  not  a  definite  point 
but  rather  a  plateau  of  high  concentration.  Figure  14  shows  the  fraction 
ef  debris  as  a  function  of  distance  from  groxmd  zero  along  the  line  of 
maximum  concentration  for  three  values  of  the  mean  vlnd  velocity.  This 
figure  indicates  that,  as  the  wind  velocity  is  Increased,  the  material  is 
spread  in  such  a  way  that  the  maximum  decreases  emd  the  small  fractions  are 
extended  to  greater  distances.  As  indices  of  the  general  natTire  of  the 
vlnd  effects,  the  generalizations  shown  in  these  three  figures  are  useful. 

They  should  not,  however,  be  accepted  as  representations  of  real  fallout 
patterns . 

If  fallout  patterns  are  to  be  forecast,  it  is  necessary  to  forecast  the 
wind  values  to  be  used  in  the  tremsport  equation.  It  is  well  known  that  wind 
forecasting  is  an  uncertain  art.  In  order  to  try  to  extrapolate  the.  degree 
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Fig.  12  —  Moximum  fraction  of  device  os  o  function 
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txC  success  which  might  be  realized  in  this  art  in  the  next  few  years,  a 
(study  of  the  application  of  numerical  forecast  techniques  for  the  fallout 
problem  has  been  attempted  at  RAND.  This  study,  Graphical  Methods  For  The 
Quantitative  Prediction  of  Close-In  Fallout,  indicates  that  useful 
bimierlcal  forecasts  can  be  made  in  some  Instances.  The  chief  difficulty, 
Which  is  peculiar  to  fallout  problems,  is  the  scale  of  motion  which  must  be 
considered.  One  of  the  chief  conclusions  in  Ref.  20  la  that  much 
more  detailed  observation  and  conputation  is  required  to  provide  the 
necessary  forecast  precision.  A  concomltteuat  of  the  small  scale  which  must 
fee  considered  is  the  removal  of  the  geostrophic  assunption  from  the  basic 
prediction  eqmtions.  A  suitable  alternative  assvmption  was  proposed  and 
tested  with  remarkable  success. 

The  forecast  scheme  used  in  Ref.  20  is  not  Intended  to  be  a  routine 
method  for  fallout  forecasting.  It  shows,  however,  that  the  problem  of 
forecasting  winds  for  the  prediction  of  fallout  may  be  soluble  by  the  proper 
application  of  modern  forecast  methods.  For  extremely  good  forecast  of  the 
fallout  patterns,  the  numericeil  approach  offers  the  best  and  most  reliable 
method  for  making  such  predictions. 
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VI.  SYNTHESIS  OF  FALLOUT  PATTERM 

The  previous  four  sections  have  taken  in  some  detail,  the  foior  major 
facets  or  ii^uts  of  the  fallout  computation  model.  This  section  describes 
the  methods  used  to  cmblne  these  basic  ccm^onents  into  a  confutation  of  the 
fraction  of  the  radioactivity  at  a  set  of  points  on  the  ground. 

It  has  already  been  assumed  that  the  distribution  of  activity  with 
height,  particle  size,  and  cloud  radius  ecre  Independent  of  one  another.  Let 
the  activity  fxinction  be  defined  so  that  the  integral  of  the  product  of  the  dis¬ 
tribution  fvinctions  over  the  entire  volume  of  the  cloud  and  over  all  particle 
sizes  is  eqTial  to  one.  Then,  if  this  function  is  integrated  over  only  those 
particle  sizes  and  Initial  heights  which  land  at  a  particular  point,  the 
result  will  be  the  fraction  per  unit  area  at  that  point.  The  transport 
equation  can  be  used  to  define  the  limits  of  particle  size  and  height  at  the 
point,  and  in  this  manner  the  problem  can  be  solved. 

Uhfortunately  none  of  the  functions  are  Immediately  amenable  to  a  sliifle 
mathematical  solution.  In  order  to  carry  out  the  indicated  integration,  two 
slnflifying  procedures  have  been  developed:  (l)  The  integration  cam  be  done 
numerically  by  defining  the  functions  ecf irically,  and  (2)  The  functions  can 
be  approximated  by  tractable  equations  and  a  mathematical  integration  can  be 
meuie.  Both  of  these  techniques  have  been  used;  the  first  is  reported  in 
Ref.  18  and  the  second  in  Refs.  19  and  21.  The  method  of  Ref.  l8  will  be 
briefly  reviewed  and  cast  into  graphical  terms  in  order  to  demonstrate  the 
procedure . 

The  first  step  is  to  confute  the  displacement  for  a  series  of  particle 
sizes,  each  from  a  number  of  different  heights,  and  to  make  a  plot  of  the 
ground  positions  of  these  particles  as  in  Fig.  I5.  Then  some  point  where  it 
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l8  desired  to  confute  the  fraction  of  the  device  Is  chosen;  call  this  point 
P(X,y).  Consider  now  a  particle  from  height  of  size  r^,  which  was  ini¬ 
tially  in  the  center  of  the  cloud.  This  will  land  at  a  distance  from  the 
selected  point  P.  Therefore,  with  the  assumption  of  no  relative  diffusion 
of  particles,  a  particle  of  size  r^^,  from  height  which  was  Initially  at 
a  distance  R^^  from  the  center  of  the  cloud  will  land  at  the  point 
The  activity  density  contained  on  particles  of  size  r^,  from  height  at  a 
distance  from  the  center  of  the  cloud  is  given  by  the  product  of  the  three 
activity  density  distributions.  If  the  activity  density  functions  for  an 
particles,  r^,,  from  all  heights,  Hj,  and  radii,  which  land  at  the  point, 
determined,  -they  can  be  integrated  numerically  over  particle  size 
and  height.  The  resulting  number  is  a  density,  or  fraction  of  device  per 
unit  area,  at  the  point  on  the  ground. 

This  Integration  forms  the  basis  of  the  computing  procedure  presented 
in  Ref.  18.  The  computational  procedure  has  been  used  to  test  the  distri— 
bxitlons  against  the  observations  of  Operations  Castle  and  Redwing  with  con¬ 
siderable  success  to  distances  up  to  about  100  mi  downwind. 

The  machine  confutation  is  ideally  suited  to  testing  the  variation  of 
parameters  for  comparison  with  observations  because  of  the  detailed  msuiner 
in  which  the  winds  are  treated.  In  order  to  investigate  the  changes  of  the 
pattern  with  wind,  however,  it  is  more  convenient  to  simplify  the  problem  so 
as  to  permit  mathematical  integration  which  will  bring  the  main  effects  of 
wind  velocity  into  clear  focus.  The  sinplif ications  invoked  were  the  use 
of  an  effective  radiological  radius  for  the  cloud,  the  assumption  of  a 
constant  concentration  with  height  over  a  layer,  and  the  approximation  of  the 
log-normal  activity— particle  size  distribution  with  an  exponential  function 
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for  particle  sizes  greater  than  25ia.  By  applying  the  mean  value  theorem  to 
the  transport  equation^  the  entire  vind  structure  i^  replaced  by  a  mean  wind. 
The  variation  of  some  salient  features  of  the  fallout  pattern  with  the  mean 
vind  could  thereby  be  coniputed  vlth  less  time  and  effort  than  that  required  by 
the  con^lex  machine  operation.  The  details  of  the  confutation  are  reported 
in  Ref.  19.  Figure  l6  is  a  graphical  representation  of  the  mathematical 
integration.  It  shows  the  fraction  of  the  device  per  unit  area  as  a  function 
Of  wind  velocity  and  distance  from  ground  zero  along  the  axis. 

In  order  to  provide  a  simplified  method  of  computation  which  will  permit 
more  detail  than  a  single  Integrated  wind,  the  basic  single  wind  model  was 
modified  to  use  integrated  winds  to  a  series  of  heights.  This  is  acconplished 
by  treating  different  layers  of  the  cloud  as  separate  entities.  The  basic 
eqviatlons  were  solved  for  six  10-kft  layers  and  the  results  of  each  calcular- 
tion  were  presented  in  graphical  form.  Scaling  laws  were  developed  to  permit 
the  graphs  to  be  used  over  a  range  of  yields  for  approximately  one  to  20  MT. 

By  making  a  plot  of  the  down  position  of  a  few  particles  (see  Fig.  15)  and 
applying  the  single  wind  method  to  each  layer  separately,  it  is  possible  to 
confute  the  fraction  of  the  device  at  a  point  by  summing  the  results  of  each 
layer.  Reference  21  presents  the  details  of  this  analysis.  Figure  17  shows 
the  results  of  applying  this  layer  method  to  the  Zuni  shot  of  Operation 
Redwing.  The  results  may  be  considered  as  an  excellent  approximation  to  the 
estimates  made  from  observations,  particularly  when  one  compares  Fig.  17 
with  Fig.  3- 
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Fig.  16  —  Fraction  of  a  3.5  MT  devices  X  10®  per  n  mi^  along  the 
axis  of  the  pattern  as  computed  from  the  single  wind  model 
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Observed  fraction  per  square  n  mi  X  10® 


Fig.  17  —  Calculated  vs  observed  fallout  for  the  Zuni 
shot  of  operation  Redwing.  Calculations  made  by 
the  six  —  layer,  hand -computing  model 
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The  foregoing  brief  svunmsiry  of  HAND'S  studies,  in  the  field  of  radio- 
active  fallout  should  point  several  salient  facts.  Perhaps  the  most 
ii^ortant  of  these  facts  is  the  absence  of  accurate  measurements  of  the 
basic  parameters  necessary  to  predict  or  compute  fallout  patterns .  The 
direct  measurements  made  to  date  do  not  provide  either  the  distribution  of 
activity  with  particle  size  or  the  distribution  of  activity  in  space.  As  yet 
there  has  been  no  atteiispt  made  to  determine  the  average  falling  velocity  of 
actual  fallout  particles.  The  restricted  number  of  test  sites  has  practically 
prohibited  any  study  of  the  effect  of  soil  type  on  the  resulting  fallout 
particles.  And,  finally,  the  current  methods  of  wind  measurement  and  fore¬ 
casting  are  inadequate  for  a  precise  attack  on  the  problem,  even  if  the  other 
factors  were  well  known. 

Despite  these  obvious  deficiencies  it  hsis  been  possible  to  construct  a 
Useful  working  model  for  the  confutation  of  the  fallout  resulting  from 
Burface-^jurst  weapons.  Further,  within  the  framework  of  this  model,  it  has 
been  possible  to  assess  the  relative  effect  of  the  varloue  parameters  on 
the  resulting  fallout  patterns.  The  most  pronounced  effect  is  the  varlatioh 
in  intensity  cavised  by  the  wind.  A  change  in  meem  wind  velocity  from  20  to 
4o  kn  cem  reduce  the  value  of  the  peak  concentration  by  a  factor  of  two  and 
increeuse  the  intensity  at  great  distances  frcm  ground  zero  by  as  much  as 
30  per  cent  (see  Fig.  l4). 

Another  parameter  which  has  a  major  effect  on  the  fallout  pattern  is  the 
height  of  the  cloud.  The  height  of  the  peak  of  the  Initial  space  distribu¬ 
tion  is  related  almost  linearly  to  the  dlstemce  to  which  the  mAvTimmi  is 
displaced  downwind  from  groxind  zero,  and  this  varies  with  yield  and,  to  some 
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extent,  vith  meteorological  conditions. 

The  distrlhutlon  of  activity  with  particle  size  is  Iniiortant,  but  not 
quite  so  critical  as  the  wind  and  the  height  distribution.  With  the  assuoqr* 
tlon  that  the  fall  rates  are  closely  approximated  by  the  conqiutatlons  for 
sphericed.  particles  of  density  2.^,  the  distribution  must  have  a  shape  some-' 
thing  like  a  log-normal  with  a  mode  approximately  at  kOn  radius,  pmall 
Variations  from  this  distribution  will  not  change  the  results  of  a  cooputtt- 
tion  markedly,  but  according  to  Pig.  2  it  would  be  i]iq>ossible  to  prodxice 
realistic  results  if  the  mode  were  cheuoged  by  more  than  lOji  either  way.  It 
is  necessary  to  note  that  the  small-particle  end  of  the  distribution  is 
Conpletely  undetermined.  Our  approach  limits  us  to  data  from  particles 
Which  fall  in  a  relatively  short  time;  thus,  the  small-particle  end  of  the 
distribution  is  pure  extrapolation. 

The  model,  together  with  possible  advances  in  the  state  of  the  wind 
forecasting  art,  offers  some  promise  of  a  reasonably  accurate  prediction  of 
fallout.  The  major  limitations  on  the  possibility  of  fallout  forecasting 
are  the  lack  of  detailed  knowledge  of  the  initial  distributions  and  our 
inability  to  predict  the  small  sc£l.e  motions  of  the  atmosphere.  These  un¬ 
certainties  shoTild  not  preclude  the  use  of  forecasting  as  a  guide  to  more 
detailed  delineation  of  the  patterns,  but  it  is  necessary  to  recognize  the 
sources  of  error  and  to  be  aware  of  the  uncertainties  inherent  in  any  fore¬ 
cast  fallout  pattern. 

Other  organizations  ha,ve  devoted  much  effort  to  the  forecasting  of  fall¬ 
out  and  have  had  moderate  success  in  forecasting  the  patterns  to  be  expected 
from  the  test  shots.  Most  of  the  methods  used  are  similar  in  principle  to 
the  methods  outlined  here,  and  the  degree  of  success  which  they  have  attained 
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in  operational  eituatlons  points  out  both  the  usefulness  of  the  fedlout 
forecast  and  its  Inherent  limitations. 

It  is  our  belief  that  the  models  developed  at  RAIH)  during  the  last  tvo 
/ears  represent  methods  of  fallout  prediction  which  are  commensurate  in  pre^ 
cision  with  the  demonstrated  ability  to  mesusure  the  pattern  of  fallout.  Ve 
feel  that  it  is  possible  to  make  io^rovemente  in  the  pure  mechanical  cooq^u- 
ting  processes  which  %d.ll  permit  operationally  iiseful  fallout  forecasts  to  be 
made.  But  we  feel  constrained  to  point  out  the  Inherent  errors  of  this  type 
of  forecast  and  to  warn  against  blind  dependence  on  any  forecast. 
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Commander  2 
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Maxwell  Air  Force  Base,  Alabama 
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Medicine 

Randolph  AFB,  Texas 
Att;  Research  Secretariat 

Commander,  1009th  Sp  Wpns  Squadron  1 

Hq. ,  USAF 

Washington  25,  D,C. 

Commander  3 

Wright  Air  Development  Center 
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Director  of  Development  Planning  2 

Office,  DCS/Development 
Hq.,  USAF 

Washington  25,  D.C. 

Att:  RAND  Project  Office 

Deputy  Chief  of  Staff,  Operations  1 
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Washington  25,  D.C. 

Office  of  Director  of  Defense  Research  1 

and  Engineering 
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Washington  25,  D.C. 
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Washington  25 j  D.  C. 

Commandant  1 
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Commandant  1 
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Chief,  Defense  Atomic  Support  Agency  15 

Washington  25,  D.C. 
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Commander,  Field  Command  1 
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Defense  Atomic  Support  Agency 
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Defense  Atomic  Support  Agency 
Sandla  Base 
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Att:  Mrs.  Bertha  R.  Allen,  Librarian 
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Operations  analysts  and  radiological -weapons -effects  people  will  be 
Interested  in  this  concluding  report  of  a  study  sponsored  by  the  Armed 
Forces  Special  Weapons  Project,  now  the  Defense  Atomic  Support  Agency.  The 
investigation  considers  the  physical  processes  of  local  radioactive  fallout 
caused  by  surface -burst  nuclear  weapons  and  the  development  of  mathematical 
models  that  would  permit  the  conputation  or  prediction  of  surface  fallout  patterns. 
(Research  memorandums  RM-1932,  RM-2006,  RM-2108,  RM-2115,  RM-2148,  RM-2193,  and 
1^-2296,  and  Classified  Paper  S-62  are  part  of  this  series.  The  basic  mechanisms 
for  understanding  local  fallout  are  also  investigated  in  R-309,  Close-In  ) 

In  liartlcular,  this  summary  discusses  the  basic  parameters  needed  to  construct 
a  fallout  computation  model.  The  parameters  examined  are  the  distribution  of 
activity  with  particle  size,  the  distribution  of  the  activity  In  space  at  the 
time  of  stabilization,  the  fall  velocity  of  the  particles,  and  the  wind  struc¬ 
ture.  The  study  considers  the  relative  Importance,  availability,  and  accuracy 
of  the  Inputs  to  the  model,  describes  the  various  RAND  fallout  models,  and 
evaluates  fallout  models  ge^ierally  as  operationally  useful  predlcition  devices. 

It  is  concluded  that  reasohably  accvurate  fallout  forecasts  are  possible  through 
the  use  of  mathematical  models,  but  that  the  operational  consumer  must  recognize 
the  sources  of  error  and  the  inherent  uncertainties  of  any  predicted  fallout 
pattern. 
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•  Report  remains  SRD,  with  Distribution  Statement  C  assigned. 
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Simulation  Techniques  for  Composite  Shield  Testing,  R.  2,  Short  Pulse  Width  (1  to  2 
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